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O ' Abstract 

I I QED Compton scattering at HERA is discussed in terms of the information it may 

^ ■ reveal on the proton structure at low momentum transfers . Detailed Monte Carlo studies 

'"^ I are performed which show that the analysis of inelastic QED Compton events allows the 

OO . extension of present HERA structure function measurements to a kinematic regime, which 

I up to now was only accessed in fixed target data. For these studies an improved version of 

i the COMPTON event generator is used, where special emphasis has been put on modelling 

\^ I the hadronic final state at low invariant masses. 

T-H ■ As the low regime is sometimes also discussed in the context of the collinear approxi- 

' mation and the possibility of measuring the photonic content of the proton, the Monte Carlo 

i studies are also used to check the validity of this approach. It is found that the proposed 

I concept of a photon density 7 of the proton docs not provide sufficient accuracy for the 

. description of inelastic QED Compton scattering. 

1 Introduction 

D ; 

• • , Measurements of deep-inelastic lepton-proton scattering (DIS) provide information that is cru- 

. ^ \ cial to our understanding of proton structure. Since the fixed target experiments have discovered 

^ ■ scaling violations [TIP] , much progress has been made in extending the kinematic regime covered 

. in terms of the Bjorken variable x and the four- momentum transfer Q^. This holds especially 

for the HERA ep scattering experiments which, with their wealth of data, have shown that 
the evolution of the proton structure function i<2(x,(5^) is well described by perturbative 
Quantum Chromodynamics (pQCD) throughout a wide range in x 

andQ2 However, 

at small deviations from pQCD predictions are observed [HlIIl, indicating the transition 
into a regime where non-perturbative effects dominate and the data can only be described by 
phenomenological models such as those derived from the Regge approach 0. 
In order to study this non-perturbative regime, the structure function F2 has been measured at 
very low values of and x, which are accessible at HERA via special devices mounted close to 
the outgoing electron beam direction |7i thus facilitating measurements of the scattered electron 
at very low angles. These devices, however, do not cover the transition region at Q'^ ~ 1 GeV^, 
which up to now has only been investigated using "shifted vertex" data [SlEj. In this paper 
the possibility to extend the kinematic domain of HERA into this region using QED Compton 
(QEDC) events, i.e. ep events with wide angle hard photon radiation, is discussed. 
The present studies are based on a modified version of the COMPTON event genera- 
tor J2j with a complete description of the low mass hadronic final states based on the SOPHIA 



*Deutsches Elektronen-Synchrotron (DESY), D-22607 Hamburg, Germany. E-mail: victor@mail.desy.de 
^Institut fiir Physik, Universitat Dortmund, D-44221 Dortmund, Germany 



1 



p 



I 




q = l - I'-k 



X 



k 



V 



P 




I 



k 
V 



Figure 1: Lowest order Feynman diagrams for the radiative process ep — > e^X with photon 
emission from the electron hne. I, P represent the four-momenta of the incoming electron and 
the incoming proton, while I', k and X are the momenta of the scattered electron, the radiated 
photon and the hadronic final state, respectively, s and t denote the squared four-momenta of 
the virtual lepton. 

model ^31) as described in section [21 Apart from studying the potential of QEDC events to 
access the low region this new version of the COMPTON program thus also allows inves- 
tigating the possibility to measure in the fixed target region at higher x, because the result of 
such a measurement crucially depends on the accurate description of the hadronic final state at 
low masses W. 

The possibility to measure the proton structure function F2 at low using QEDC events 
was first discussed by Bliimlein et al. |14|ll5j. In the framework of the equivalent photon 
approximation the authors introduced the concept of a photon density of the proton^, 7, to 
be valid at very low virtualities Q^. This function has been computed by de Riijula and 
Vogelsang when proposing an extraction method for 7 from HERA QEDC data. The 
validity of this approach is discussed in the second part of the paper, where it is shown that the 
collinear approximation does not provide a sufficient description of inelastic QEDC scattering 
and that measurements of the QEDC cross section with reasonable precision can thus not be 
interpreted in terms of the photon density function. 



Radiative processes in ep scattering, as depicted in Fig.^ may be split into three different 
classes |17(ll8j with (i) the bremsstrahlung or Bethe-Heitler process corresponding to small 
masses of both the virtual electron and the virtual photon, (ii) the QED Compton process 
with a low virtual photon and a large virtual electron mass and finally (iii) the radiative DIS 
process where the photon is collinear either with the incoming (Initial State Radiation, ISR) 
or the outgoing (Final State Radiation, FSR) electron. All three classes correspond to distinct 
experimental signatures. For the QEDC scattering process the final state topology is given by 
an azimuthal back-to-back configuration of the outgoing electron and photon detected under 
rather large scattering angles. In this configuration their transverse momenta balance such that 
very low values of the exchanged photon virtuality are experimentally accessible. 

To correctly describe the process ep e'yX the standard kinematic variables x and Q^, used 
to describe inclusive deep-inelastic scattering (DIS), have to be redefined in order to account 

^Sometimes also denoted as /^|p or D^^p. 



2 QEDC Monte Carlo Simulation 
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Figure 2: Energy of the photon candidate hi the backward calorimeter as calculated by the 
HI simulation program. Shown are the expectations from QEDC scattering as predicted by 
the COMPTON generator in comparison to the DIS background simulated using DJANGO. 
The left and right plots show the distributions before and after applying the cut of 9(f on the 
maximum polar angle of any further energy deposition in the calorimeters (see text). 
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Here I and P are the four-momenta of the incoming electron and the incoming proton, and I' 
and k represent the momenta of the scattered electron and the radiated photon, respectively 
(Fig.^. Three further independent variables are needed for a full description of the differential 
QEDC scattering cross section. In the formalism presented in the Lorentz invariant scale 
variable x-y = q ■ l/P ■ I and the scattering solid angle defined in the centre-of-mass frame 
of the virtual Compton process and encapsulating two degrees of freedom are employed. The 
cross section is then given by [T7] 



dxdx^dQ'^dO.* 



/-)/* /pi-^y ■^7j Q 



da 

In* 



da 
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where [da/d^*]'^'^ are the differential cross sections of the process ej* — > for transverse 
and longitudinal polarised photons, fully calculable in the framework of QED ^7j, and /^'y^ 
represent the corresponding virtual photon spectra, which may be expressed in terms of the 
photo-absorption cross sections cj^*p. In order to specify a'J^ip one has to consider three separate 
contributions depending on the value of the invariant mass W of the outgoing hadronic final 
state: 



1. Elastic scattering, for which the proton stays intact (W = mp). This channel is well 
measured, and the cross section is given by the electric and magnetic form factors Ge 
and Gm', 

2. Resonance production, where the total mass of the hadronic final state X lies in the range 
mp + m^ <W <2GeY; 

3. Continuum inelastic scattering at > 2 GeV. In this region the 7*p cross section is 
defined through the proton structure functions F2 and F^. 

The above cross section expression (j2| is implemented in the COMPTON event generator |12j . 
However, as this generator was primarily written for an application in analyses of elastic QEDC 
events, in the original version a rather crude approach has been employed to describe the 
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Table 1: Summary of QEDC selection criteria as described in text. 



Item 



Cut value 



QEDC signature 



E^,E^ >4GeV 

Ee + E^> 20 GeV 

A = 180° - A(j)< 45° 

E > 0.5 GeV (within one cluster) 



Hadronic final 
state 



9 < 90° 



Event properties 



existence of reconstructed vertex 



resonance region and only simple scale invariant F2 parameterisations are used to model the 
continuum inelastic domain. Furthermore, no hadronisation of the final state X is performed. 
As this paper aims at the investigation of inelastic QEDC events a new version of the COMP- 
TON generator was developed which includes detailed parameterisations for the reso- 
nance and the continuum pUj regions. In addition, several packages for a complete simula- 
tion of the hadronic final state have been implemented into the program. For the present studies 
the SOPHIA package ^Hl is used in the range of low and low masses, W, of the hadronic 
final state while the Quark Parton Model with subsequent Lund string fragmentation 21^ is 
employed at high W and high Q^. 

3 On the Measurement of F2 Using QEDC Scattering 

The presented studies are based on a Monte Carlo sample which was generated using the 
new version of the COMPTON program and corresponds to an integrated luminosity of about 
30 pb~^; the incident beam energies used are E^ = 27.6 GeV for the electron and Ep = 820 GeV 
for the proton beam. In order to obtain a realistic simulation of the experimental conditions at 
the HERA ep detectors, the generated events were subject to the GEANT-based simulation of 
the HI detector Hence, the selection criteria are adapted to the resolution and acceptance 
limits of the detector components relevant for this analysis. The most important ones are the 
backward calorimeter^ SpaCal [231 ? liquid-argon (LAr) calorimeter and the tracking and 
vertex detectors [22| l2l|. 

The analysis strategy is based on the detection of the outgoing photon and the outgoing electron 
observed almost back-to-back in azimuth. Thus, the main experimental requirement is the exis- 
tence of two electromagnetic energy depositions (clusters) with E^, E^. > 4 GeV in the backward 
calorimeter; an additional limit is imposed on the total 67 energy Ee + E^ > 20 GeV to reduce 
radiative corrections. In order to separate QED Compton scattering from the FSR process and 
still obtain inelastic QEDC events the azimuthal acoplanarity angle A = 180° — Ac/) is required 
to be below 45°, where A4> represents the angle between the transverse momenta of the electron 
and the photon; this cut also matches a corresponding requirement in the COMPTON event 
generator J7|- Inelastic QEDC events are selected demanding at least one LAr cluster with 
energy above E > 0.5 GeV; as the acceptance of the LAr calorimeter is limited to ^ > 4° this 
cut — apart from rejecting elastic events — substantially reduces the inelastic contribution at 
very low masses W. After applying these criteria, background contributions from standard DIS 
events only appear for large angles of the hadronic final state particles where one hadron fakes 
a photon signal in the backward calorimeter and the current jet is observed in the backward 

■^The z axis of the right-handed coordinate system used by HI is defined to he along the direction of the 
incident proton beam and the origin to be at the nominal ep interaction vertex; the backward direction is thus 
defined through 2; < 0. 
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Figure 3: Kinematic domain 
of continuum inelastic QED 
Compton events in compari- 
son to the regions covered by 
inclusive DIS measurements 
at HERA and fixed target ex- 
periments. 
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region of the detector. Tiiis contribution is reduced by requiring no additional energy deposi- 
tion apart from the two clusters to be found above ^max = 90° (see Fig.|21). Finally, only events 
with a reconstructed vertex are considered, in order to guarantee a correct determination of 
all kinematic variables. However, as the final state hadrons are generally scattered under very 
small angles the corresponding tracks do not provide sufficient resolution for an accurate vertex 
measurement. Vertices are thus reconstructed using the electron track detected in one of the 
tracking detectors. A brief summary of all selection criteria is given in Tab. ^ 
The phase space covered by a QEDC sample selected by the cuts described is shown in Fig.|31for 
an integrated luminosity of 30 pb~^. Compared to the kinematic range accessed at HERA via 
standard deep- inelastic scattering, the QEDC events clearly extend to lower Q^. For inclusive 
DIS the outgoing electron is not detected for such low values of as it is scattered at small 
angles escaping through the beam pipe unseen; QEDC events, however, with the electron and 
photon in the final state balancing in transverse momentum, reach into the transition region 
below < 1.5 GeV^, which otherwise is only accessed through data taken in the shifted vertex 
runs piirU] or with dedicated devices for tagging events at very low like the ZEUS BPT Q. 
But, due to acceptance limitations in the high x domain, these latter data cannot extend the 
low F2-iiieasurements to such high x as QEDC events. It is therefore the range of medium 
to high X which is of special interest when analysing QEDC scattering. As higher x correspond 
to low masses W, special emphasis had to be put into the correct modelling of the hadronic 
final state. 

An accurate description of the hadronic final state is especially important as, for the kinematic 
range in question, the reconstruction of the Bj or ken- variable x cannot be performed using the 
kinematics of the outgoing electron and photon; the x-resolution of this method deteriorates 
with 1/y = xs/Q"^ thus becoming inapplicable at low values of the inelasticity y. For a double 
differential measurement of the structure function F2{x,Q^) at 1 GeV and x > 10~^ 

the variable x has thus to be reconstructed from the final state hadrons. This is done using 
the so-called Sigma method [25;, which is based on the measurement of ^i{E — pz)i summing 
over all objects i of the hadronic final state. To suppress the infiuence of calorimeter noise 
substantially contributing at low masses W, we use a simple approach in which only clusters 
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Figure 4: Correlation between the generated and reconstructed variables and x using the 
Sigma method as described in the text. 



with energies above the noise level of 0.5 GeV are considered when reconstructing x. Fig.l^shows 
the correlation between the generated and reconstructed values of and x. It demonstrates the 
possibility to reconstruct the QEDC event kinematics using the hadronic final state throughout 
the relevant phase space region, i.e. for 0.1 < < 10 GeV^ and 10^^ ^ ^ ^ 10^^. 
The expected statistical significance of an F2 measurement based on the same COMPTON 
event sample with 30 pb~^ total integrated luminosity is presented in Fig.|SJ In order to extract 
the structure function F2 in terms of x and the selected Monte Carlo data are divided into 
subsamples corresponding to a grid in x and Q^. The bin sizes are adapted to the resolution in 
the measured kinematic quantities such that purity and stability in all bins shown are greater 
than 30%; here, the purity (stability) is defined as the ratio of the number of simulated events 
originating from and reconstructed in a specific bin to the number of reconstructed (generated) 
events in the same bin. As the proton structure function F2 is then obtainted by a bin-by-bin 
unfolding method using the same Monte Carlo sample, the extracted F2 values trivially lie on 
the curve representing the ALLM97 F2 parameterisation used as input to the Monte Carlo 
generator. Clearly, the systematic uncertainties remain to be studied in order to judge on the 
overall precision of the measurement. 

The results demonstrate that an analysis of QEDC scattering data at HERA would add informa- 
tion in the low and medium to high x region not yet covered by the HERA ep experiments. 
Such an analysis would extend the present HERA measurements into the region previously 
covered only by fixed target data. 



4 On the Photon Content of the Proton 



In contrast to the exact treatment of the QEDC scattering process, the concept of the pho- 
ton content in the proton based on the collinear or equivalent photon (Weizacker- Williams) 
approximation ^niEOl provides a much simpler approach to QEDC scattering and is believed 
to reveal basic features of photon- induced reactions involving proton beams. In this "parton 
model" approach the transverse component of the exchanged photon momentum is neglected 
and the emitted photon is assumed to be on-shell and collinear with the incident proton. This 
simplifies the expression for the QEDC cross section to jH^ 



2 



where the "structure" function 'y{xi, Qf) parameterises the photon-parton content of the proton 
depending only on two degrees of freedom, the leptonic variables Qf = 1 — V and xi = 2pf^ j 
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Figure 5: Possible measurement as expected from a 30 pb^^ QEDC event sample. The 
expected HERA results (closed circles), estimated using the COMPTON Monte Carlo pro- 
gram, are shown in comparison to HI results obtained from standard DIS data (open trian- 
gles l^) taken during the years 1996-1997 and measurements from fixed target experiments 
(open squares |26y . open stars \27\ and open crosses .28]). 



fixed centre-of-mass energy ^/s the inelasticity yi is also defined through these two variables via 
the relation yi = Qf /xis. 

According to the above expression, an experimental determination of the double differential 
QED Compton scattering cross section d?a/dxi dQf can be interpreted as a measurement of 
the photon-parton density function ^{xi,Qf). This 7 function can then be applied for the 
computation of other ep and pp cross sections where the underlying process is mediated by 
a quasi-real photon exchange [HSIISSIIS^IHSIISSIISZI- The possibility to measure the structure 
function ^(xi,Qf) using QEDC events was discussed in several publications [ni[T3l[T6..38) . We 
consider here the latest work by De Rujula and Vogelsang JHl where they proposed to compare 
the Qf dependence of 7 at fixed x/ with its gluon counterpart g(x,Q'^). 

In order to judge on the feasibility of measuring 7 a dedicated study of the accuracy of the 
collinear approximation is performed. Fig.^lshows the transverse momentum {pt,-y*) distribution 
of the exchanged photon in elastic and inelastic QEDC events as predicted by the COMPTON 
event generator. While for elastic events pt^-y* is rather small, much larger values are reached 
when selecting inelastic scattering processes. Thus, as in the collinear approximation one as- 
sumes to vanish, one can expect significant deviations when calculating event kinematics. 
The effect is further enhanced due to the acceptance constraints, which demand that both out- 
going particles, electron and photon, are measured under finite polar angles. For the presented 
studies the QEDC selection criteria listed in Tableware applied on generated quantities; in addi- 
tion the detector acceptance is approximately reproduced by demanding 0.06 < 9e, 0^ < tt— 0.06. 
These criteria are close to those proposed in (IB) . 
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Figure 6: Transverse momentum of the exchanged photon in elastic (left) and inelastic (right) 
QEDC events generated hy the COMPTON program. 



A comparison between the exact calculation and the predictions from the collinear approxima- 
tion is given in Fig. [3 which shows the total, i.e. the sum of the elastic and inelastic QEDC 
cross section as a function of Qf in bins of x/. Significant discrepancies are observed in several 
bins. The exact cross sections were derived from a sample of COMPTON events generated 
without radiative correction. For the collinear approximation predictions have been computed 
analytically by Vogelsang "39^ using the same set of cuts; the corresponding uncertainty of this 
calculation was estimated to be approximately 2%. Furthermore the elastic scattering cross 
sections generated by the COMPTON program and computed in the Weizsacker- Williams ap- 
proximation were compared. Very good agreement was obtained showing that the collinear 
approximation provides, indeed, a good description of the elastic process. 

It has been checked that the difference between the F2 parameterisations employed in the 
COMPTON generator and in the calculations by De Riijula and Vogelsang, cannot account for 
the observed discrepancies: the different predictions for 7 calculated according to the relations 
given in jl5ll6j agree within 4% ^. The inelastic QEDC cross sections, however, deviate by up to 
a factor of 2 from each other. This shows, that apparently the equivalent photon approximation 
does not provide a sufficient degree of accuracy in the inelastic region and that the factorisation 
of the cross section given by Eq. ^ in terms of only two kinematic variables is a too rough 
approximation of the QEDC scattering process. 

In order to ensure the validity of the collinear approximation, additional selection cuts were 
proposed in jTH] : 

-i,s >lGeV'^ and pt,e,Pt,^ > IGeV , (4) 

where the momentum scales t and s represent the virtualities of the exchanged lepton in the 
Feynman diagrams shown in Fig.^ Here, the cuts applied on the transverse momenta pt^e and 
Pt^ry, actually, have a much stronger effect on the distributions of COMPTON events than those 
imposed on t and s. After applying all of the four additional cuts the COMPTON QEDC cross 
section is again compared to the corresponding analytical calculations provided by Vogelsang. 
As shown in Fig. |H1 there before, significant discrepancies observed. 

Fig. ISl illustrates the difficulties of the collinear approximation. For Qf < 5 GeV^ and 5.6 TO^^ < 
XI < 1.8 • 10"^ the cross sections predictions from the COMPTON program drop significantly 
when applying the additional cuts on t, s and pt^e, Pt,f This is not the case for the equivalent 

•^As the phenomenological F2 parameterisations used in the COMPTON program have, contrary to the pQCD 
based PDFs employed in I16| . no artificial lower limit, we also tested different start scales Qq for the integration 
over F2, as specified in |lfci) and Eq. (1) of |15|. The differences in 7 are in any case of the order of a few percent. 
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photon approximation where the cross sections remain unchanged due to an incorrect calculation 
of the event kinematics arising from the asumption pt^-y* = 0. 

It is not ruled out that some carefully chosen cuts on the transverse momenta, polar angles and 
other event quantities may limit the phase space such that the equivalent photon approximation 
becomes applicable also for inelastic QEDC scattering. However, with such limitations it is 
questionable whether the 7-function describes a characteristic and universal property of the 
proton relevant for precision measurements. 

5 Concluding Remarks 

The presented analysis shows that QEDC scattering data at HERA provide the potential to 
measure the proton structure function F2 at low and medium to high x, extending the 
kinematic range covered so far into the fixed target regime. As such a measurement requires a 
good understanding of the hadronic final state at low invariant masses, an improved version of 
the COMPTON generator was developed which allows for a better event modelling. With this 
generator the accuracy of calculations using the equivalent photon approximation is studied. 
It is revealed that this approximation is not able to describe the inelastic QEDC cross section 
with sufficient precision. A measurement of the photon density of the proton from QEDC events 
would, thus, not provide any decisive insight into the proton structure. 
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Figure 7: Double differential cross section of the QED Compton scattering. Open circles depict 
cross section values given by the COMPTON generator. The corresponding error bars show 
statistical errors. The dashed lines denote the values computed by W. Vogelsang. 
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Figure 8: QED Compton scattering cross section after additional kinematic cuts. Open circles 
depict cross section values given by the COMPTON generator. The corresponding error bars 
show statistical errors. The dashed lines show the values computed by W. Vogelsang. 



12 



